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ABSTRACT 

Using canonical (Schrodinder) quantization of spherically symmetric gravi- 
tational dust systems, we find the quasi-classical (coherent) state, 
that corresponds to the classical Schwarzschild solution. We calculate the 
"quasi-classical Schwarzschild metric" , which is the expectation value of the 
quantized metric in this quasi-classical state, gffi =< a^'lg^a^ >. De- 
pending on the quantization scheme that we use, we study three different 
quasi-classical geometries, for all of which gj%) turns out to be singularity 
free. Their maximal extensions are complete manifolds with no singularities, 
describing a tower of asymptotically flat universes connected through Planck 
size wormholes. 
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Spherically symmetric dust matter in 4D Einstein gravity is an inte- 
grable system [1]. One can fix the (diffeomorphism) gauge completely, and 
get an infinite set of independent (non propagating) degrees of freedom [2]. 
Though very simple, this dust system is rich enough to have some phys- 
ically interesting (classical) solutions. Two of these are the homogeneous 
dust-filled Robertson- Walker universe and the Oppenheimer-Snyder collaps- 
ing star. The first gives us a simple model of the big bang, and the second 
is a simple model of black hole formation. Because of the simplicity of this 
dust system, it can be easily quantized [2,3]. In this letter we study some of 
the effects of this quantization. Especially, we find the quasi-classical metric, 
and study its structure. 

Following Lund [2], we first fix the gauge classically, and then quantized 
the reduced system. One can fix the gauge by choosing comoving coordinates 
in which the metric has the form 



where \i and A are functions of r and 71, and dQ| is the volume element in 
S 2 . Explicitly, we choose r and 1Z such that 



where U u are the 4-velocities of the dust particles, and the primes denote 
"^.-differentiation. This gauge fixing corresponds to N± — 1 , N l — , and it 
is complete [1,2]. The resultant reduced Hamiltonian [2,3] is 



where p = y/l — 1Z 2 , p s is the surface of the dust star, Y(t, p) is the (reduced) 
field variable, Y = e A , and P Y = 5L/5Y = 3ttYY /(Zp^/T^p 1 ). In the 
(r,p,9,(f>) coordinates, the metric (1) is 



ds 2 = _ dr 2 + e 2, dn 2 + e 2X dn 2 



(1) 



(2) 




(3) 



ds 2 = -dr 2 + (d p Y) 



2 dp 2 
1-p 



2 + Y 2 dVi\ 



(4) 
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The Hamiltonian (3) is time independent, namely Ti{p) = E(p), where 
E{p) is a function of p only. Classically, E(p) is determined by the Cauchy 
data, Y(t — 0) and Y(r — 0). A most important feature of (3) is that there 
are no p-derivatives in it. This means that the infinite number of degrees 
of freedom are decoupled. Each can be described as an harmonic oscillator 
[3]. To see this it is more convenient to use a different gauge. Consider the 
gauge N± = Y/p and N p = — J v d p (Y/ p)dr( , which corresponds to choosing 
the time coordinate r\ instead of r, where 



J p 



(5) 



The metric in this gauge is 



ds 2 



Y 2 

-drf 
p 1 



{d P Yf T 



d p (Y/p)d V > 



dp 2 



d p (Y/p) 2 d V ' 



dr]dp + Y 2 dtt 2 2 



(6) 



where (d p ) r denotes p derivative in which r is held fixed, while d p (with no 
t index) denotes a p derivative where rj is fixed. Redefinition of the field 



X(r,,p) = Y(r,,p) 



3tt 



1 - ?E(p) 



(7) 



implies that the following modified Hamiltonian, H, is 77-independent 

rps y/l - p* 



H 



\ 6ir 



p2 



3tt 



X z dp 



E 2 (p)dp (8) 



where Px = 3^(8^1 — p 2 ) 1 {d v X) 2 . Clearly, the system described by (8) is 
an infinite set of (decoupled) harmonic oscillators^. 

2 To see that, one can discretize p by denning pk — jfPs, k = 1, 2, N. Eq. (8) is then 
replaced by a set of TV (one dimensional) harmonic oscillators, satisfying 

H k = auPl^) + k Xl(r}) = E k 



where X k (r]) = X(rj,p k ), P k (rj) = P(ij,p k ), E k = E(p k ), and a k 
In the continuum limit, N — > 00, this is an infinite set. 



6tt ' P k ~ 



:i7r 
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Considering classical solutions, we choose r = (or equivalently 77 = 0) 
such thatfl P(rj = 0, p) = 0. In this case, the general classical solution of (8) 
(or (3)) is 

X(7 1 ,p)=X (p)cos(r ] ) (9) 

where X (p) = X{jj = 0,p). For a given initial data, X (p), or equivalently 
for a given initial dust density pdust(v = 0, /?), we have a unique solution (9). 

There are two interesting (classical) solutions: the homogeneous density 
solution, and the zero density one. For the homogeneous density case, we 
have 

X = xt om \p) = la p and Y^ hom \r, p) = l -a p{\ + cost?) (10) 
where is a constant and r = |ao(?7 + sin (77)). The metric (4) is now 

dt? = _ dr 2 + /I + CQS ^\\jf + p 2 d ^ ) (n) 



2 v w/ 7 l-p 

which is the well known cycloidal solution for a homogeneous dust-filled 
universe. 

For the zero density solution we have 

X^X { t h \p) = j 2 (12) 

or 

y(^)(r,p) = ^(l + cosfa)) 
P 

M 

t = -^(v + sin(v)) (13) 

Here M is a constant, and we will see that it is indeed the Schwarzschild 
mass. The uniqueness theorem, tells us that the only spherically symmetric 
zero density solution is the Schwarzschild one. To see that (13) is indeed the 
Schwarzschild solution, we write it in Schwarzschild coordinates (t,r,9,4>). 

3 Of course this is possible only for a certain class of solutions, but we are interested 
only in this class. 
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The metric in the (r, p) coordinate is (4) with (13), from which we see that 
the Schwarzschild radial coordinate, r, should be defined such that 

r = Y {Sch \r 1 ,p) (14) 

Define the Schwarzschild time, t, such that g# = —l/g rr , to get the familiar 
form 

(d.,%„, = -(l~) * 2 + Y^mJ;. + r'dal (15) 

An interesting feature of these two solutions, is that they can be matched 
continuously at the surface of the dust star, p = p s [4,5]. From (10) and (13) 
we see that the matching condition isQ 

2M = a Q pl (16) 

This is the Oppenheimer-Snyder (OS) collapsing star [4]. 

The coordinates (r, p, 8, <p) are the Novikov coordinates, where our p is 
related to the original Novikov coordinate, R*, by p 2 = (R* + 1)~ 2 . The 
range of p is < p < p s < l.In the OS solution, p runs twice over the interval 
[0, p s ]; from to p s to cover the interior homogeneous region, and then back 
from p s to 0, to cover the exterior Schwarzschild region^. 

The quantization of this simple dust system is elementary. It is easily 
done using (3) or (8). In coordinate representation 

X(p, v ) = X(p, V ) , P x ( p , v ) = -th^ ( i7) 

we get the quantum version of Eq. (8), which is the Schrodinger equation 

= j E(p)<S>\X; n ]dp (18) 

4 Together with (14), that says r s (r) = a(r)p s . 

a In the outside (Schwarzschild) region, p runs only over the semi-open interval (0,p s ]. 
It does not reach p = which is asymptotic infinity r = oo. 
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where m(p) = — — and u> = 1. We use a Schrodinger quantization 

scheme, and not a Wheeler-DeWitt one. Namely, the scalar product of quan- 
tum states (which determines the Hilbert space) is the Schrodinger oneQ, 
< ^i|^2 >= / ipl(x)ipz(x)dx. The solutions of (18) are harmonic oscillators, 
and the spectrum E(p) is quantized [3]. In the homogeneous case this leads 
to a quantization condition for clq [6,8,9], and in the collapsing OS case, it 
leads to a discrete mass [9,10]. 

The metric (6) is now an operator (replace Y with the operator Y). We 
would like to calculate the expectation value of this quantized metric in a 
quasi-classical (coherent) state, |a[X(p);^] >. We call the resultant matrix 
of real functions 

9$(v>P)=<<x\9i*>\<x> ( 19 ) 

the "quasi-classical metric", which corresponds to the classical solution (9). 
Using (18), we find the quasi-classical states, which are the eigenstates of the 
annihilation operator 

a x \a[X{p)- V ] >=a[X(p);r ] }\a[X(p);r ] } > (20) 

where 

and Xq(p) is the classical initial value of X. As in the standard harmonic 
oscillator case, the coherent state, \a >, is not and eigenstate of the Hamil- 
tonian (18). It is a superposition of such states. But one can give the usual 
statistical interpretation, it describes the probability of finding the system 
(collapsing star or R-W universe) with a given mass (or energy). The prob- 
ability is described by a very sharp Gaussian around the classical value. 

6 The Wheeler-DeWitt scalar product is of a Klein-Gordon type [6,7] 
< *i|* 2 >= / 4>*( x ) d x ifa(x) . 
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Knowing the quasi-classical states (20), and the metric operator which in 
the Heisenberg representation is 



Y(n,p) = X (p)l + X(n,p) 



—tr) 



(21) 



we can calculate the quasi-classical metric. One should notice however that 
we have an ambiguity in the quantization scheme that we can use. The 
commutation relations of a x (p) and aJ x (p) can be functions of p. And different 
functions correspond to different quantization schemes. This is related to the 
ordering problem in our Schrodinger equation [6,7]. At this point we don't 
need to determine the explicit ordering that we use. We define the following 
creation and annihilation operators 



a = 



m(p)f(p) 



-a(p) and a) 



I 



m(p)f(p) 



d\p) 



(22) 



where the commutation relations of a and a* are p-independent, [a, at] = 
1. The p-dependent commutation relations of a x (p) and a x (p) is given by 
m(p)f 2 (p). So the function f(p) determine the quantization scheme (the 
ordering) that we use. With this definition we get the following metric ex- 
pectation value 



< aids \a > 




dn 2 + 



- d. 




rj 



dp 2 - 









dndp + fdnj } (23) 



P z 



where (ds 2 ) is (6) with the classical Y(n,p) solution (9). As expected, the 
quantum corrections are proportional to H, and depend on the quantization 
scheme (f(p))- 

Consider first the homogeneous case. There is only one degree of freedom 
in that case [2,9] , it is the scale factor, a(n). In that case Aa and AP a should 



7 



notbe functions of p. This corresponds to choosing /(p) 
p is a constantQ In that case we get [1 1] 



pp in (22), where 



<a {h) \ds z \a w > 



CO 



M + ^ + P 2 dnl\ (24) 



This describes a ball of dust that collapses to a minimum radius (the Planck 
radius) and re-expands again [11]. The quasi-classical metric (24) is singu- 
larity free. 

In the Schwarzschild case the situation is more complicated but as we are 
going to see, one can avoid the classical singularity also in that case. Using 
the Schwarzschild radial coordinate (14), we see that as long as r 2 >> hf 2 , 
the quasi-classical metric (23) is almost the classical Schwarzschild one. The 
deviations are of order hf 2 /r 2 or lpf 2 /r 2 . As we expected, the quantum 
deviations are very small for r > lp. On the other hand near the classical 
singularity, r = 0, the quantum corrections are dominant. From (23) we have 



K) 



Sch 



(2Mr . , 2 ^ (d pV ) 2 T 



+ o(r 2 ) 



l-p 2 



+ 



+ 



{(d P f) 2 + (d pV ) 2 T f) ^ ~ (d P (f/p))W + o(r 2 ) 
" (^(d P (f/p) 2 )v + o(r 2 )) d V dp + fy 2 + o(r 2 )) dtf 2 , 



dp 2 
(25) 



near r = 0. We see that the dominant parts are the quantum corrections. To 
study the structure of the metric (25) (near r = 0) we need to know (d p r)) r . 
One should remember that t(i], p) is not the classical function (12), there are 
quantum corrections. To find r we write the metric (23) in the form 



ds„ = —dr 



g rp drdp + g pp dp 2 + ( Y 2 lass + \f 2 



dn 2 2 



(26) 



and solve for r(r/, p). Generally it is quite complicated to get r explicitly, it 
is a complicated function of the classical solution and /(p). But near r = 



7 The uncertainty principle is (Aa)(AP a ) = \h, so p 2 > 1. 
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one can easily solve for r and get 




+ r 



(27) 



and 



i]d p log 



f{p) 



P 



) 



(28) 



as r 



0. 



In the classical Schwarzschild case {d p 7]) T diverge like r _1 as r goes to 
zero. This is because of the classical singularity. In the quasi-classical case 
we see (from (28)) that generally (d p r]) T is finite as r goes to zero. From 
the uncertainty principle we know that f(p) cannot be zero (in the inter- 
val (0, ps\), so (28) cannot diverge for a C^-function, f(p). If f(p) is also 
bounded, the Riemann tensor will be everywhere finite, and generally its 
maximum will be around r = and of order Mjl\. So for a general quan- 
tization scheme (for which f(p) is a bounded C 1 function), one avoids the 
classical Schwarzschild singularity. Because there is no singularity at r = 0, 
one can extend the space-time manifold to negative r's. Depending on the 
function f(p), we may have different situations: From (13) and (14) we see 
that r = corresponds to r\ = ir, so the metric (25) on r = (and fix (#, </>)), 
is 



In this paper we examine only the cases in which r[/] does not change sings 
in the interval (0,p s ], so we have three different situations: i) r[/] > 0, in 
which case r = is spacelike, ii) r[/] = 0, in which case r = is null, 
and iii) r[/] < 0, for which r = is timelike. Typical Penrose diagrams of 
the maximally extended manifolds for cases (i), (ii) and (iii) are shown in 
Figures 1, 2 and 3, respectively. In Fig. 1 the manifold contains regions (7) 




(29) 



where 




(30) 
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r > r h , (II) < r < r h , (II*) -r* h < r < 0, and (I*) r < -r h . Where 
r h = 2M + o(M 2 /M 2 ), and -r* = -2M + o(M 2 P /M 2 ). Regions J* and II* 
are very similar to I and JJ respectively, but of course r — > — r (as well as 
i — ► — i) there. In Fig. 2 the manifold contains the regions (I) r > rh, (II) 
< r < r h and (///) r < 0. In this case, r = is a Cauchy horizon like 
the inner horizon of a Kerr black hole. In Fig. 3 the manifold contains the 
regions (/) r > r + = rh, (II) r_ < r < r + , (III) < r < r_ and (IV) 
r < 0. In this case, since r = is timelike there is also an inner horizon ( 
at < r_ < r + = rh ) in which ds 2 a = 0. This manifold is very similar to a 
Kerr one, but with no singularity at r = (or anywhere else). 

The manifolds in Figures 1, 2 and 3 are all complete. Namely, in all 
of which all timelike curves have infinite proper length. They all describe 
asymptotically flat universes (r — ► ±oo) which are connected through Planck 
size "wormholes" (r = 0). 

From (14) we see that r is the expectation value of the dynamical field 
operator, r =< a\Y\a >. In a standard quantization one should consider 
both possitive and negative field configurations, so the quantum extension 
to negative r's seems to be "natural". 

In this letter we use the midisuperspace of spherically symmetric dust 
universes, suggested by Lund, to study the quasi-classical structure of the 
quantized Oppenheimer-Snyder collapsing star. Especially, we find the out- 
side quasi-classical Schwarzschild metric. Depending on the quantization 
scheme that we use, the classical singular hypersurface r = 0, becomes a 
regular spacelike, null or timelike hypersurface. The quasi-classical metric is 
generally everywhere singularity free. The maximal analytic extension of this 
space-time is a complete manifold, with no singularities. It describes a tower 
of asymptotically flat universes connected through Planck size wormholes. 
The cases in which r = is not always spacelike, null or timelike[|, require 
further investigations. 

Can we go beyond this midisuperspace model? Clearly, when we include 
non-spherical deviations, we get propagating graviton degrees of freedom, 
and the quantization become problematic. It could be though, that non per- 
turbative approaches may give consistent results [12], otherwise we need to 

8 This corresponds to sings changes in T[f] (see eq. (29)). 
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modify Einstein gravity or quantum physics. On the other hand, it may be 
interesting to extend Lund's midisuperspace but still consider only spher- 
ically symmetric geometries. For example, spherically symmetric charged 
dust systems. 

There is however one general feature that one can draw. Our dynamical 
variable, Y(i],p) is exactly the vierbein, e^(x). In a general theory for which 
the vierbein is the basic dynamical variable, the metric operator is its square, 
9fiu = e^e aj ,. If we use a free field representation for the vierbein field, we 
can describe it as a collection of harmonic oscillators, and have the standard 
quasi-classical coherent state. The metric expectation value in this quasi- 
classical state will be < g^ v >= (e^e au ) c ias+o(h) , and the classical singularity, 
= 0, can be avoided. 
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Figure Captions 

Fig. 1: Perose diagram for a typical maximally extended manifold in the 
case where r = is spacelike. This is an infinite tower of maximally extended 
Schwarzschild spacetimes. Each two successive Schwarzschild manifold, M 
and M*, are related by r — > — r and t — > —t. 

Fig. 2: The Penrose diagram of the maximally extended manifold in the 
case where r = is null. The classical singularity r = 0, becomes a Cauchy 
horizon. 

Fig. 3: Penrose diagram for a typical maximally extended manifold in the 
case where r = is timelike. The manifold is very similar to the maximally 
extended Kerr manifold, but there is no singularity at r = 0. 
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Figure 1: 

Penrose diagram for the maximally extended manifold where r=0 is spacelike. 
This is an infinite tower of maximally extended Schwarzschild space-times. 
Each two successive Schwarzschild manifolds ( M and M*) are related by 
r --> -r and t --> -t. 




Figure 2: 

Penrose diagram for the maximally extended manifold where r=0 is null. 
The classical singularity, r=0, becomes a Cauchy horizon. 




Figure 3: 

Penrose diagram for the maximally extended manifold where r=0 is timelike. 
The manifold is very similar to the maximally extended Kerr manifold, but 
there is no singularity at r=0. 



